Observation of wave turbulence in vibrating plates 



Arezki Boudaoud 
Laboratoire de Physique Statistique, CNRS / ENS / Paris 6 / Paris 7, 
24 rue Lhomond, 75231 Paris Cedex 5, France 



00 

o 

O 
(N 

> 

o 
O 



Q 
U 



> ■ 
00 . 

O ; 

d 

00 

o 



Olivier Cadot, Benoit Odille, and Cyril Touze 
ENSTA-UME, Unite de Recherche en Mecamque, 
Chemin de la Huniere, 91761 Palaiseau, Cedex, France 
(Dated: October 6, 2008) 

The nonlinear interaction of waves in a driven medium may lead to wave turbulence, a state such 
that energy is transferred from large to small lengthscales. Here, wave turbulence is observed in 
experiments on a vibrating plate. The frequency power spectra of the normal velocity of the plate 
may be rescaled on a single curve, with power-law behaviors that are incompatible with the weak 
turbulence theory of Diiring et al. [Phys. Rev. Lett. 97, 025503 (2006)]. Alternative scenarios are 
suggested to account for this discrepancy — in particular the occurrence of wave breaking at high 
frequencies. Finally, the statistics of velocity increments do not display an intermittent behavior. 

PACS numbers: 47.27.Gs, 62.30.+d, 47.35. Jk 



The statistical distribution of energy and energy fluxes 
are central questions concerning out-of-equilibrium dissi- 
pative systems with a large number of degrees of free- 
dom. When waves propagate in a medium, their nonlin- 
ear interaction might generate other waves with different 
wavenumbers, which means that energy is transferred be- 
tween different lengthscales. If the amplitude of waves 
is large enough, this transfer leads to a distribution of 
energy on a large number of wavelengths, and the sys- 
tem reaches a state called wave turbulence [T] , such that 
the energy cascades between scales and might be dissi- 
pated on a small scale. Although they share the same 
phenomenology, wave turbulence is much more advanced 
analytically [1] than hydrodynamic turbulence 2]. For 
waves of small amplitude, the framework of weak tur- 
bulence yields kinetic equations, the solutions of which 
have been derived starting from the mid-1960s and cor- 
respond to energy spectra with power-law dependence on 
the wavenumber. Wave turbulence might apply to capil- 
lary ^] or gravity waves on the surface of liquids, 
to plasmas 31j to nonlinear optics to magnetohydro- 
dynamics [Sj] or even to Bose-Einstein condensates llCf]. 

Experimental studies are much less numerous than 
theoretical ones; they were performed either on the 
oceanographic scale — waves on a stormy sea (e.g. [lll|). 
or on the laboratory scale — capillary and gravity 



waves fll [H, Q El in, [13, [3- Besides, the domain 
of validity of weak turbulence theory is still a matter of 
debate. On the one hand, discontinuities in the slope 
of breaking waves result mathematically in a wide en- 
ergy spectrum 1^, 23], as apparently observed for grav- 
ity waves [l3]- On the other hand, weak turbulence the- 
ory results in Gaussian statistics for the waves, in con- 
trast with experiments when bursts of intense motion 
occur 13|, llTI ] , a phenomenon known as intermittency. In 



as it provides a new system, vibrating plates, where wave 
turbulence could be observed. 

Here we study experimentally a suspended plate driven 
at high amplitudes [ll]. We show that a wide energy 
spectrum is generated, discuss its interpretation in terms 
of weak turbulence and wave breaking, and investigate 
whether the system is intermittent. The typical broad- 
band spectrum observed is also of special interest for its 
applications, e.g. for reproducing the sound of thunder 
in theaters. It is also related to the bright shimmer- 
ing sound of gongs and cymbals [2^, [3] ■ Transition to 
chaotic vibration was studied for cymbals in [23|] , and for 



panels in [25|, |2 



this context, the theoretical study in [2l| is very useful 



The experimental setup consists of a steel plate sus- 
pended to a rigid frame and forced with a vibration gener- 
ator (shaker B&K4810, glued to the plate with beeswax) 
moving perpendicularly to the plate (fig.[T^). The plate 
comes from a reverberation unit named EMT140, that 
was widely used in studio recordings to add a reverber- 
ated sound effect to dry signals recorded by near-field 
microphones 27[. Hence, the plate was chosen for its 
very high modal density, obtained by large dimensions 
2 m X 1 m for a thickness oi h — 0.5 mm, as well as 
for the moderate values of the quality factor, in order to 
get a fuzzy reverberated sound. Material properties were 
estimated as: Young's modulus E = 200 GPa, Poisson's 
ratio V = 0.3 and mass per unit volume p = 7800 kg/m'^. 
The plate is fixed at its four corners, so that the bound- 
ary condition is mainly free. The forcing is sinusoidal at 
/i = 20 Hz that is close to a resonant frequency of the 
plate; it was chosen in order to allow the best injection 
of energy in the system, so that the turbulent regime is 
reached more easily. A laser vibrometer gives the normal 
velocity v(t) at a given point in the plate. The signal is 
acquired at the sampling frequency of 32 kHz, and the 
FFT is computed from 50 s of signal, averaged over time 
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FIG. 1: The experiment, (a): set-up with a steel plate of 
dimensions 2 m x 1 m and thickness h — 0.5 mm; close-up 
view of the fixation. (b),(c): time series of the local transverse 
velocity measurements v{t) for the forcing frequency fi — 20 
Ifz; duration of 10 times (b) and 1 time (c) the forcing period. 



windows of 0.5 s, so that A/ = 2 Hz. A force sensor 
(impedance head B&K 8001) is mounted between the 
shaker and the plate. The simultaneous measurement of 
the velocity at the same point gives the average power 
/ = {Fv) injected by the generator into the system (with 
1 mW of accuracy) . 

For a bending wave of frequency / and wavenumber fc, 
the dispersion relation is 

/ = hck^, with c = ^E/12p{l - z/2)/ 2tt (1) 

proportional to the sound velocity in the bulk material. 
It was checked in ^27| that this dispersion relation indeed 
holds in the present setup. It gives the space-time cor- 
respondence of the statistical properties of the velocity 
signal, similarly to Taylor's hypothesis for fully turbu- 
lent flows [2], when fluctuations are not too large. 

For very low forcing amplitude, the velocity signal v{t) 
recorded by the vibrometer is sinusoidal. For higher am- 
plitude, it becomes chaotic (fig. [T]d,c). In the frequency 
space, v{t) is characterized by its power spectrum Pv{f), 
given by the Fourier transform of the auto-correlation 
function, P^{f) = J {v{t)v{t + t)) exp{2TrifT) dr. This 
spectrum becomes broadband at high forcing (fig. [2la,), 
which is typical of wave turbulence; however, even with 
a long time-averaging of the signal, keeps a num- 
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FIG. 2: Power spectra of the transverse velocity, (a): raw 
Pv{f) as a function of frequency /, for different values of in- 
jected power I (in increasing order as displayed by the arrow: 
< 1 mW, 2.3, 8.8, 26.4, 68.8, 136 mW); errors A/ = 2Hz and 
APv = 10~^ (mm/s)^/Hz. (b): rescaled spectra {Io/I)Pv ac- 
cording to Eq[2]vs. f/fc for all forcing amplitudes, where fc 
is defined by {Io/lf^^Pv{fc) = 10"^ mm^s. Inset (c): evo- 
lution of fc with the forcing intensity . The continuous line 
is the best power law given in Eq. 0, yielding an exponent 
a = 0.33. 



ber of peaks corresponding to the plate eigenfrequen- 
cies. We checked homogeneity (by changing the exci- 
tation and measurement points) and independence on 
boundary conditions (by imposing fixed displacements at 
points at the edge); these changes affected very slightly 
the power spectra below the injection frequency . 

As the forcing amplitude is increased, the spectra ex- 
hibit a wider and wider power-law dependence on fre- 
quency, Pv{f) ^ / with /? = 0.5 ± 0.2 (this error is an 
upper bound), which would correspond to the cascade 
regime. It is followed by a fall which could correspond to 
the dissipative scale. We seek the best rescaling of the 
spectra as a function of the injected power /. This yields 
the scaling form 

F.(/) = (///o)l/^0(///c), /c CX /,;(///o)". (2) 

Here a scaling function, Iq a unit of power and fc a 
cut-off frequency. This rescaling enables to collapse the 
spectra on a single curve (FiG.[2h)- The exponent for the 
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FIG. 3: Damping factor for the present plate, from [27 



dependence of /c on / is found to be a = 0.33 ± 0.01 
(fig. [2t). In the cascading frequency range, this imphes 

Pv{j) f-P = /0-66±0.07j-0.5±0.2 

In order to compare with previous theoretical work, we 
first note that the power spectrum for the transverse dis- 
placement ^ of the plate is given by P^{f) oc P-u(/)A£?. 
When weak turbulence is attained, as investigated in [2]| , 
the spatial power spectrum of the displacement can be 
rewritten as P^{k) oc c^^e^^'^k~'^ , introducing energy 
flux e per unit mass (e has units of a velocity cubed 
and is proportional to the power input /), and omit- 
ting numerical prefactors and a logarithmic dependence 
on k. e is proportional to the injected energy /. The 
1/3 exponents for e comes from the ^ — s- — ^ symme- 
try of the plate, which involves four waves interaction. 
The spectrum can be translated into the frequency space 
P({k)kdk oc P^(/)d/. Using the dispersion relation ([1]), 
we get P^{f) cx he^/^f~'^ and Pv{f) oc he^/^ is constant. 
This dependence is significantly weaker than in the mea- 
surements (FIG.[2^,b). 

In the framework of weak turbulence, nonlinearities 
of order p imply that Pvif) scales as e^/^ {H; the ex- 
ponent 1/2 obtained for p = 2 is the closest to the 
measured .66 ± 0.07 ([2]). This value of p = 2 means 
three waves interactions, a quadratic nonlinearity and 
no ^ — > — ^ symmetry. Indeed, geometrical imperfec- 
tions are unavoidable in real plates, which is known 
to break this symmetry and to produce quadratic non- 
linearities [2^. Therefore, we assume in the following 
that Pvif) ~ e^^^/~^^^7 corresponding to a displacement 
spectrum P^{f) ^ e^/^f-^f^. This assumption allows to 
investigate the possible role of damping in setting the 
cutoff frequency. 

Indeed, we introduce the damping rate j{f) ^ ■ 
The spectrum of the energy per unit mass is £{f) oc 
Pvif) ^ e^/^/^^/^. Let us consider the balance of 
energy over the cascade frequency range; the in-flux 
is e while the energy dissipated till fc — hckc is 
f^-fik)£{k)dk oc /^=7(/)f(/)\/W7d/. Balancing 
these two fluxes yields fc ~ e^^"^^. For our setup, a 



flt to the damping coefficient measured in [27[ is shown 
in FIG. [3] in the frequency domain of interest. It yields 
6 ~ 1/2, so that fc e ^ I , which is far from measure- 
ments ([2|) and so damping cannot account for the cutoff. 

A last option is that the wide energy spectrum might 
be generated by singularities of the plate displacement as 



for gravity waves [l9|. For plates, wave breaking would 
be replaced 21 1 by ridges 2^, 30| and d-cones [sil, [l^l ■ 
It was shown in [20| that random independent slope dis- 
continuities result in a spectrum P^{f) oc i/sF^/"*, I's 
being the frequency of occurrence of slope discontinuities 
and r the rms velocity impulse at each discontinuity. For 
the velocity Pvif) oc z/^F^/"^ which compares with the 
second part of the spectra (fig. [2k,b) over half a decade. 
Besides, the jump should be given by the typical rms 
velocity Urms so that we expect F ~ e^/^. As a con- 
sequence, the whole spectrum could result from a 3- waves 
interaction for low frequencies, as suggested above, and 
singularities for higher frequencies. These two spectra 
match at a frequency f = fc such that ei/^/"!/^ 
yielding / e^/"^, which agrees with the scaling ^ as 
seen in fig. 

Finally, we consider the statistics of of the velocity in- 
crements defined as ArV = i;(i + T) — i;(i). The PDFs are 
displayed in fig. for the large forcing amplitude. An 
intermittent behavior of the velocity statistics would be 
revealed by a change in the PDFs shape as the lag r de- 
creases . Here we can see in fig. that the PDF shape 
remains satisfactorily Gaussian whatever r. The struc- 
ture functions, Spir) = (|Ai-w|p), are plotted in fig. [4)d. 
They are generally used to determine the scaling behavior 
of the velocity differences statistics with the time-lag r 
0. The structure functions start to decrease for r < 50 
ms (i.e. the forcing period). For very small r < 0.3 
ms (i.e. the cut-off frequency), the velocity signal be- 
comes smooth and a simple scaling behavior Spir) = 
is found. For wave turbulence, the range of interest is 
comprised between these two last extremes. However, 
within this range no clear power laws are distinguishable 
in FIG. [4)3. We then chose to plot the structure functions 
versus 5*2 (t) in fig. Hj;. This technique was used for fully 
developed turbulence to measure anomalous scaling ex- 
ponent due to the intermittency phenomenon [33| . In our 
case, the scaling exponents, defined as: Spir) oc S'2(t)'»p, 
are indicated in fig. [H; for each order moment p. There is 
no significant deviation from (^p = p/2, meaning that no 
anomalous scaling is observable. Hence, wave turbulence 
in plates does not exhibit any intermittency phenomenon. 

To summarize, we observed a broadband spectrum 
in a vibrating plate and investigated the variations of 
the cut-off frequency. In this context, internal damping 
mechanisms (mainly thermoelastic and viscoelastic losses 
for our plate [13]) seem to be irrelevant. Losses at the 
edge [21J can be discarded as the plate is fixed only at 
the corners. The radiation of acoustic waves in air is 
negligible since the frequencies of interest are well below 
the coincident frequency, for which bending and acous- 
tic waves have the same phase velocity. The value of 
this frequency has been measured as 20 kHz in our set- 
up [27] . For thicker plates, the coincident frequency may 
fall in the frequency range of interest, thus leading to a 
huge increase of the damping factor, see 6.17. [34| . This 
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FIG. 4: Statistical properties of the velocity increments 
AtIi = v(t + t) — (injected power 136 mW). Probability 
density functions compared to Gaussians in (a). Structure 
functions of order p = 2,3,4,5,6, versus (b): the timelag r, 
and (c): the order 2 structure function S2. Continuous lines 
are best power laws fits with exponents C,-p (see text). 



could affect the conclusions on the cut-off. Our experi- 
mental results suggest a 3-waves spectrum matched to a 
spectrum of singularities where dissipation occurs. Ob- 
viously they call for more theoretical effort, in particular 
concerning the weak turbulence of plates with quadratic 
nonlinearities or the turbulence of singularities. 
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